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I.  INTRODUCTION 


A.  MOTIVATION 

Close  proximity  towing  of  surface  ships  is  a  matter  of  interest  to  the  U.S.  Navy 
and  the  Office  of  Naval  Research.  Several  possible  applications  have  been  suggested 
including  the  SLICE/KAIMALINO  connection  (Nash)  and  the  SEA  LANCE 
configuration  (TSSE).  Advantages  of  ship  towing  are  cost  effectiveness  when  towing 
large  expensive  loads  with  a  denser  vehicle  and  being  able  to  tow  various  kinds  of 
platforms  to  carry  various  kinds  of  loads  (Nash).  One  of  the  problems  that  need  to  be 
overcome  is  the  evaluation  of  motions  of  the  two  ships  when  in  close  proximity  towing 
operations.  This  thesis  addresses  this  problem,  and  suggests  a  possible  design 
methodology. 


B.  PREVIOUS  WORK 

Nash  utilized  strip  theory  calculations  (Beck  and  Troesch)  in  order  to  develop  the 
coupled  heave/pitch  equations  of  motion  of  the  two  ships.  A  brief  discussion  on  strip 
theory  is  available  in  “Principles  of  Naval  Architecture,  Volume  III,  Chapter  VIII, 
Section  3.4.”  The  matching  condition  that  was  used  was  that  of  a  generic  spring/damper 
connector.  In  this  thesis,  we  want  to  utilize  a  more  realistic  matching  condition  based  on 
the  speed/resistance  profile  of  the  ships.  In  order  to  avoid  the  excessive  computational 
burden  of  fully  coupled  hydrodynamic  and  structural  analyses  (Korsmeyer  and  Kring), 
we  wish  to  develop  a  simpler  design  procedure  that  can  be  effectively  utilized  in  design. 
This  procedure  could  be  applied  in  cases  where  only  the  operational  envelopes  of  the 
ships  are  known,  and  not  necessarily  their  individual  hydrodynamic  properties. 
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II.  REGULAR  WAVES  RESPONSE 


Beck  and  Reed  describe  ship  seakeeping  as  one  of  the  most  challenging  problems 
in  the  marine  hydrodynamics  field.  The  complexities  of  wave  resistance,  maneuvering 
problems  and  unsteadiness  due  to  incident  waves  as  well  as  some  of  the  nonlinearities 
associated  with  free  surface  problems,  coupled  with  the  need  for  easy  implementation  in 
a  design  environment,  contribute  to  the  challenges  we  face  in  seakeeping  assessment. 

In  our  study,  we  assumed  that  the  water  is  incompressible,  the  density  is  constant 
and  the  flow  is  irrotational.  These  are  standard  assumptions  in  seakeeping  studies.  We 
also  assumed  that  the  ships  operate  in  infinitely  deep  water. 

In  this  chapter,  the  waves  that  the  ships  are  encountering  are  assumed  to  be 
regular  waves.  This  motion  of  waves  is  two  dimensional,  sinusoidal  in  time  with  radian 
frequency  co .  The  wave  propagates  with  velocity  V  such  that,  to  an  observer  moving 
with  this  velocity,  the  waves  appears  as  steady-state  (Newman).  Wave  surface  elevation 
can  then  be  defined  as 

r/{x,  t)  =  cos{kx  -  co  t)  (1) 

where  the  parameter  k  =  co/V , which  is  called  the  wave  number,  the  number  of  waves 
per  unit  distance  along  the  x-axis.  Clearly,  k  =  27rl  A,  where  the  wavelength  A  is  the 
distance  between  two  points  on  the  wave  with  the  same  phase.  For  infinitely  deep  water, 
the  dispersion  relation  states  that  k  =  co2  /  g  (Newman). 

The  frequency  of  the  wave  encounter  by  the  ship  is  defined  as 

coe  =a>  +  Vk  cos  /?  (2) 

where  j3  is  the  relative  wave  heading  and  V  becomes  the  ship  speed  relative  to  water 
(McCreight). 

We  can  derive  the  equations  of  motion  for  free  oscillations  of  a  rigid  body  moving 
in  sea.  Motion  of  that  body  can  be  described  in  six  degrees  of  freedom  system.  Three 
translational  (surge,  sway,  and  sway),  and  three  rotational  (roll,  pitch,  and  yaw)  motions 
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are  required  for  this  system.  Figure  1  illustrates  these  motions  and  the  standard  sign 
convention  used  in  seakeeping  calculations. 


Vi 


171  =  surge  1)3  =  heave  1)5  -  pitch 

in  =  sway  i)\  —  roll  ije  =  yaw 


Figure  1 .  Sign  Convention  for  Ship  Motions  (From:  Beck) 

From  Newton’s  second  law  of  motion,  a  translational  displacement  is  defined  as 

Afj  =  F  (3) 

where  A  is  the  apparent  displacement  of  the  body  and  F  is  the  total  force  acting  on  the 
body.  For  rotational  motion  we  have 

Irj  =  F  (4) 

where  I  is  the  mass  moment  of  inertia  about  that  axis  and  F  is  the  total  moment  acting  on 
the  body. 

For  the  simplified  case,  the  total  force  and  moment  consist  mainly  of  hydrostatic 
and  hydrodynamic  forces 

F(t)=Fex(t)+FH(t)  (5) 

where  Fex  is  the  exciting  force  due  to  waves  and  Fh  is  the  radiation  force  due  to  motion  of 
the  ship. 


4 


The  exciting  force  can  be  defined  in  frequency  domain  as  Fex(t)=Fexe1,Uet . 

For  sinusoidal  motions,  the  hydrodynamic  radiation  force  and  moment  can  be 
defined  as 

F^  =  -[A  +  B  (co)fj  +  C  (co)tj]  (6) 

where  A  is  the  added  mass,  B  is  the  hydrodynamic  damping,  and  C  is  the  restoring  force 
or  moment. 

Since  our  equation  is  linear,  it  is  usually  carried  out  in  the  frequency  domain  for 
easier  calculations.  Now  we  can  define  the  motion  and  its  derivatives  in  the  frequency 
domain  as 

f](t)  =  . 

/;(/)  =  -o)e2rjei">e‘ 

Like  a  two  degrees  of  freedom  system,  we  can  write  our  equation  of  motion  as 

\~co e2  (A  +  A)  +  i afi  +  C]^  =  ¥ex  (7) 

where  (A+A),  B  and  C  are  the  mass,  damping  and  stiffness  matrices  respectively. 

For  six  degrees  of  freedom  system,  this  equation  has  to  be  defined  for  six 
different  motions. 

Now,  let’s  define  the  equations  of  motions  for  a  point  on  SLICE  and  a  point  on 
KAIMALINO. 

[A,  +  As  ]  {ij}  +  [Bs  ]  {T]}  +  [CS]  {rj}  =  {. F }  +  {/,} 

[A*  +  Ak  ]  {ij}  +  [Bk  ]  {ff}  +  [CK  ]  {7}  =  {Fex  }  +  {fK}  .  (8) 

where  fs  and  f  K  are  the  connection  forces  acting  on  SLICE  and  KAIMALINO 
respectively,  where  f K  =  —fs  .  The  connection  force  and  motions  are  illustrated  in  Figure 
2. 
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SLICE  p 
hs 


Figure  2.  Connection  Force  and  Motions  of  the  Ships 

As  the  derivation  procedure  is  the  same  for  both  equations,  let’s  work  on  a 
generic  one.  When  we  define  the  equations  in  the  frequency  domain, 

H*.2  (A  +  A)  +  i  coeB  +  C]  {7}  =  {FJ  +  {/} .  (9) 

Since  the  exponential  exists  in  all  tenns,  it  is  canceled  and  the  equations  of 
motion  in  the  frequency  domain  become  as  follows.  Let’s  define 

A  =  \—(Oe  (A  +  A)  +  icoeB  +  C] . 

We  know  that  A  is  a  6x6  matrix  for  six  degrees  of  freedom  system. 

a7j  =  {f} 

rj  =  inv{A){F)  (10) 

where  {A}  =  {Fex }  +  {/} . 

Ship  motions  due  to  regular  waves  of  a  given  wavelength  and  direction  are  now 
detennined  for  a  given  forward  speed  (V).  Motions  in  vertical  and  horizontal  planes  are 
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usually  considered  as  decoupled  and  may  be  solved  as  two  distinct  3x3  systems  vice  the 
6x6  system.  We  performed  our  studies  on  the  vertical  plane  response  of  the  ships. 
Therefore,  only  heave  and  pitch  responses  and  the  resultant  interactions  between  SLICE 
and  KAIMALINO  in  these  three  degrees  of  freedom  will  be  analyzed  from  this  point 
forward.  Hence,  all  motions  except  r|3  and  r|5  are  set  to  zero.  The  expanded  equations  of 
motion  in  two  degrees  of  freedom  become  as 


heave  - 
pitch  - 


33 


53 


A 

A 


35 


Vs 

_ 

A+f 

Vs 

Fs+fxs 

(ID 


Now  we  can  define  the  equations  of  motion  in  vertical  plane  (but  in  two  degrees 
of  freedom)  for  two  vehicles  as 

As,sTJ3S  +A35,sr/5'S  =  F3S  +fs 


j4-53,sT/3s+A55,sTJ5S—F5s  fsXS  (12) 

A33,kT/3k  +A35,kJI5k  =F3k  +fK 
As3,kT/3  k  +  Ass.kMs.K  ~  F5  k  ~fsXK 

where  x  is  the  distance  between  the  connection  point  and  the  center  of  gravity  of  that 
vehicle. 

In  order  to  solve  these  four  equations,  we  make  the  following  substitution. 

Vs ,S  ~  Ms ,S  ^3 ,sfs 

'  •"  ■  ’  ■ '  (13) 

V3 X  ~  M3  ,K  +  ^3,k/k 

Vs  ,K  ~  Ms,K  “*"^5  xfn 
where  fs  =-fK  =-/. 
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Referring  to  Figure  2,  we  can  form  the  absolute  motion  at  the  connection  points  as 


M  3,5  Ms,SXS 
,K  ~Ms,KXK 

When  we  combine  (13)  and  (14)  as 

‘zs  —  M3 ,S  —  ^3,sf  —  Ms,SXS  ^ 5,SXsf 
~  Ml,K  +  V3,nf  ~  Ms,KXK  ~V5kXkI 

Let  a  in  terms  of  n  ’s  and  b  in  terms  of  V  ’s  be  defined  as 


(14) 


(15) 


a  s  M3,K  Ms ,KXK  Ms,SX 

b  =  _ */3,.S'  ~^3,K  "F  ^5,KXK  Ft'j  jX 

Therefore,  our  equation  becomes  as 


s 

s 


Is -4k  =a~bf- 


Again  referring  to  Figure  2, 
f  =  T- 


_t$s~Ak 


l 


(16) 


(17) 


(18) 


Combining  (17)  and  (18)  gives  us  our  final  equation,  which  gives  us  the 
connection  force  between  two  ships  becomes  as 


f  =  T—-  0*) 

—  +  b 
T 

The  tension  at  the  connection  is  assumed  to  be  the  tension  at  steady  state.  It 
actually  varies  with  the  motions  of  the  ships,  but  we  assume  that  the  influence  of  its 
variation  on  heave/pitch  motions  is  of  higher  order  and  can  be  neglected.  Figure  3  shows 
the  tension  curve  we  used  in  the  calculations.  This  is  a  curve  fit  of  actual  experimental 
data  for  the  trailing  ship,  the  KAIMALINO.  It  is  of  course  valid  only  for  the  range  of 
speeds  shown  in  the  Figure. 

The  conditions  of  the  regular  waves  are  chosen  to  be  varying  speed,  heading  and 
towline  length.  Typical  results  are  presented  in  Appendix  A,  Figures  A-l  through  A-15. 
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In  all  figures,  motions  are  dimensionless  with  respect  to  a  unit  wave  height,  while  forces 
are  dimensionless  with  respect  to  earth’s  gravitational  constant,  ship  length,  water 
density,  and  area,  as 


F 

dim  ensianless 


F 

PgAL2 


(20) 


where  p  =  1.9905  lb! ft3,  g  =  32.2  ft! sec2,  Z  =  105/tand  A  =  \ft2 . 


For  consistency,  the  leading  ship  (SLICE)  length  was  used  as  the  reference  length 
(105  ft)  in  the  non-dimensionalization  of  the  results.  Zero  heading  corresponds  to 
following  seas,  while  head  seas  correspond  to  180  degrees. 

The  figures  of  the  motions  of  the  ships  show  us  that  when  the  ships  are  in  tow, 
they  move  virtually  together  at  the  corresponding  connection  points.  This  becomes 
progressively  more  accurate  as  the  towline  length  is  reduced.  This  proves  the  accuracy  of 
the  codes  generated. 

Based  on  the  results  presented  in  Appendix.  A,  we  can  draw  the  following 
conclusions. 
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1 .  The  magnitude  of  the  connection  force  is  of  similar  order  as  the  magnitude  of  the 
wave  exciting  forces.  Therefore,  its  effect  cannot  be  neglected,  and  in  fact  it  will 
influence  the  resulting  motions  as  much  as  the  exciting  waves. 

2.  The  connection  force  becomes  smaller  for  larger  towline  lengths  as  expected. 

3.  Ship  motions  at  the  connection  points  are  closely  matched  for  small  towline 
lengths,  while  we  observe  some  variation  for  larger  towline  lengths.  This  is 
consistent  with  physical  intuition. 

4.  The  dependence  of  the  connection  force  on  forward  speed  is  not  monotonic.  Even 
though  there  is  a  large  variation,  no  significant  trend  is  evident.  A  similar 
conclusion  holds  for  ship  motions  in  terms  of  forward  speeds. 
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III.  RANDOM  WAVES  RESPONSE 


We  assumed  that  the  waves  were  regular  waves  in  Chapter  II.  If  an  observer  tries 
to  observe  the  waves  encountering  his,  or  her,  observation  point,  it  is  not  difficult  to  see 
that  waves  in  real  life  are  random  both  in  space  and  in  time.  Hence,  it  is  really  difficult  to 
model  ocean  wave  behavior.  In  addition  to  this,  waves  show  different  characteristics  in 
different  regions  of  the  world. 

If  our  observer  keeps  a  record  (primarily  wave  height  and  period)  of  the  waves  for 
that  particular  point,  and  sums  those  regular  waves,  a  wave  system  for  that  point  and  for 
that  particular  time  is  gained.  We  call  a  seaway's  "spectrum"  as  the  probabilistic  function 
developed  by  taking  the  Fourier  transform  of  the  correlation  function  for  free  surface 
elevation  (Nash).  The  correlation  factor  contains  the  data  our  observer  recorded.  The 
spectrum  S(oj)  is  a  measure  of  the  energy  contained  within  a  wave  system.  Figure  5 
shows  a  typical  wave  spectrum  and  the  algebraic  addition  of  the  waves.  The  area  under 

the  curve  represents  the  mean  energy  stored  in  a  particular  wave  system,  E  =  J  S{co)dco 

o 

(Cummins). 


Figure  4.  Wave  Spectrum  and  Addition  of  Waves  (From:  Cummins) 


There  are  several  spectra  that  are  used  in  random  wave  computations  to  compute 
the  ship’s  response.  In  our  case,  the  Pierson-Moskowitz  spectrum  is  chosen  as  a 

representative  spectrum.  This  two-parameter  model  predicts  the  wave  spectrum  for  fully 
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developed,  long  crested  seas  with  no  swell  from  other  generating  areas.  Fully  developed 
seas  contain  waves  at  equilibrium,  independent  of  fetch  and  duration  of  wind.  Long 
crested  seas  have  parallel  crests  and  are  assumed  to  be  unidirectional. 

The  Pierson-Moskowitz  spectrum  is  described  by 

2 

S(<y)  =  ^exp[-/?(g/Fw)4].  (21) 

CO 

where  g  is  acceleration  of  gravity  [cm/sec“],  co  is  wave  frequency  [radians/sec],  a  is 
8.10xl0'3 ,  (3  is  0.74  and  Fw,  is  wind  speed  [cm/sec]  at  19.5  m  above  the  surface. 


Another  parameter  considered  in  random  wave  calculations  is  the  significant 
wave  height  //l/3 ,  which  is  the  average  of  the  highest  one  third  of  all  the  waves 
(Newman).  Using  an  empirical  relationship  between  wind  speed  and  significant  wave 


height,  HU3  =2 


V., 


g 


a 

Jj 


the  Pierson-Moskowitz  spectrum  now  becomes  a  function 


of  significant  wave  height  as 


,  8. MO  3 

S(co)  = - - — exp 

co 


-.032 


f  \2 

g 

Hv 

\  A  J 


co 


(22) 


Also  known  as  the  motion  transfer  function,  the  response  amplitude  operator 
(RAO)  maps  the  complex  response  of  a  vessel  to  a  seaway  or  input  spectrum  as  a 
function  of  frequency.  The  spectrum  of  the  particular  response  that  the  RAO  has  been 
computed  for  is  calculated  from: 

S^^RAOicof  S{co).  (23) 

In  this  equation,  SR  ( CO )  is  the  response  of  the  vessel  to  the  input  sea  spectrum  for  a  given 
frequency.  Ship  response  then  can  be  calculated  for  a  given  wave  frequency  and 
significant  wave  height.  For  example,  the  complex  absolute  motions  predicted  in  regular 
wave  modeling  (£,s,  qO  are  converted  into  RAO’s  for  absolution  motion  as 
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RAO<lt)  =  abs(lt). 


(24) 


The  response  spectrum  for  absolute  motion  then  becomes  as 

•s«_£,,(®)=|«te(£,,)|25<®)-  <25> 

For  our  case,  the  connection  force  response  can  be  defined  as 

S„  _ S(a>)  .  (26) 


One  way  to  predict  a  measure  of  performance  over  a  range  of  frequencies  is  by 
calculating  the  significant  double  amplitude  of  the  response.  Double  amplitudes  are 
obtained  by  integrating  the  response  with  respect  to  frequency  over  the  frequency  range 
of  any  input  spectrum.  For  our  case,  the  significant  double  amplitude  of  the  absolute 


motion  of  SLICE  at  the  connection  point  is<T^  =  j  5^  _<^s(CO)dCO .  Now  we  can  make 


(On 


more  accurate  calculations  to  make  a  better  design  with  random  wave  response. 


Different  conditions  with  respect  to  speed,  heading,  towing  length  and  significant 
wave  height  are  chosen  for  the  results.  Appendix  B.,  figures  B-l  through  B-40,  shows 
typical  results  for  these  conditions.  As  in  chapter  II,  all  results  are  non-dimensionalized 
as  shown  in  equation  (20)  and  the  related  paragraph.  In  some  cases,  designated  by  f/F3 
the  connection  force  is  further  normalized  with  respect  to  the  heave  wave  exciting  force. 


The  results  are  in  root  mean  square  (RMS)  form.  This  is  very  convenient  for 
further  calculations.  Multiplying  the  RMS  value  by  two  gives  us  the  significant  wave 
amplitude,  which  is  the  average  of  the  highest  one  third  of  the  amplitudes  of  the  response. 
Also,  multiplying  the  RMS  value  by  four  gives  us  the  eight-hour  maximum  value,  which 
can  be  used  as  design  extreme  amplitude  of  the  response. 


Based  on  the  results  presented  in  Appendix  B.,  we  can  draw  the  following 
conclusions: 


1.  Connection  forces  are  higher  in  higher  sea  states  as  expected.  Head  seas  produce 
the  highest  values  of  connection  forces. 
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2.  Relatively  (with  regards  to  wave  exciting  forces)  high  values  of  the  connection 
force  exist  in  headings  other  than  head  seas.  This  may  be  significant  in  operations, 
since  even  though  the  actual  value  of  the  connection  force  may  be  “small”  it  may 
be  considerably  higher  that  the  wave  exciting  force  and  may,  therefore, 
significantly  limit  the  operability  window  of  the  ships.  This  item  is  further 
addressed  in  the  next  chapter. 

3.  The  connection  forces  are  not  a  monotonic  function  of  the  towline  length,  as 
Figure  B-16  demonstrates.  Therefore,  it  is  possible  to  select  a  towline  length  so 
that  the  connection  force  in  minimized. 

4.  In  general,  higher  forward  speeds  will  generate  higher  RMS  values  for  the 
connection  force,  although  again  the  trend  is  not  monotonic.  As  Figure  B-36 
shows,  it  is  possible  to  select  a  forward  speed  so  that  the  connection  force  is 
minimum. 


14 


IV.  SEAKEEPING  EVALUATIONS 


Pierson  and  Moscowitz’s  sea  state  table  is  used  to  obtain  the  significant  wave 
heights.  The  mean  wave  heights,  used  as  the  representative  significant  wave  heights  in 
MATLAB  computations  for  different  sea  states,  are  shown  in  Table  1. 


SEA  STATE 

MIN. WAVE 

MAX.  WAVE 

MEAN  WAVE 

HEIGHT 

HEIGHT 

HEIGHT 

1 

0.50 

1.20 

0.85 

2 

1.50 

3.00 

2.25 

3 

3.50 

5.00 

4.25 

4 

6.00 

7.50 

6.75 

5 

8.00 

12.00 

10.00 

6 

14.00 

20.00 

17.00 

7 

25.00 

40.00 

32.50 

Table  1.  Significant  Wave  Heights 


Having  obtained  all  the  necessary  responses  in  a  seaway,  we  want  to  introduce  a 
performance  degradation  factor  to  be  used  in  design.  The  question  we  have  to  ask 
ourselves  is:  If  a  ship  can  operate  in  some  condition  (speed,  heading,  sea  state)  without 
towing,  at  what  condition  can  she  operate  at  undertow?  The  approach  we  select  is  to  fix 
the  speed  and  for  each  heading,  calculate  the  “equivalent”  sea  state.  This  equivalent  sea 
state  is  defined  on  the  basis  of  heave  and  connection  exciting  forces.  We  call  two  sea 
states  as  “equivalent”  for  our  analysis,  if  they  produce  the  same  total  vertical  force  on  the 
ship.  This  will  be  the  heave  exciting  force  for  the  condition  without  towing,  and  the  sum 
of  the  heave  exciting  force  plus  the  vertical  connection  force  for  the  condition  with 
towing.  Both  values  are  in  terms  of  their  RMS.  The  fact  that  we  use  their  algebraic  sum 
designates  an  envisioned  worst-case  scenario  where  the  phase  angles  of  the  two  sets  of 
forces  are  the  same. 
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The  results  are  presented  in  Appendix  C.,  figures  C-l  through  C-56.  The  sea  state 
variation  legend  in  these  figures  indicates  the  equivalent  sea  state  (outer  envelope)  for  the 
ships  under  towing  while  the  inner  curve  is  the  corresponding  sea  state  without  towing. 
The  sea  state  degradation  legend  in  the  figures  is  the  one  that  can  be  used  more  directly  in 
design.  It  maps  a  given  sea  state  condition  (outer  envelope)  into  an  equivalent  less 
severity  sea  state  condition  under  towing  (inner  envelope).  This  new  sea  state,  together 
with  the  vertical  connection  force,  will  match  the  severity  of  the  initial  sea  state  without 
towing.  Therefore,  we  can  use  these  results  to  map  an  existing  operability  envelope  of  a 
ship  without  towing  to  an  equivalent  operability  envelope  of  the  ships  under  towing.  A 
notional  example  of  such  an  operability  envelope  is  shown  in  Figure  5.  The  triangular 
area  shown  in  the  figure  has  been  mapped  to  the  dark  area  shown.  Using  this  procedure, 
we  can  map  any  existing  operability  area  of  the  ships  between  free  running  and  towing 
conditions. 


180 


90 

8 


Figure  5.  A  Notional  Ship  Operability  Envelope  Mapping. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  main  conclusions  from  this  work  can  be  summarized  as  follows: 

1.  Connection  forces  under  tow  are  comparable  in  magnitude  to  wave  (heave) 
exciting  forces. 

2.  The  developed  sea  state  degradation  factor  may  be  a  viable  solution  to  assessing 
ship  operability  envelopes  under  tow. 

3.  The  developed  procedure  can  be  used  to  maximize  the  expected  operability 
envelope  in  a  combination  of  sea  states. 

B.  RECOMMENDATIONS 

As  recommendations  for  further  research,  we  suggest  the  following: 

1.  What  is  the  effect  of  tension  variations  on  vertical  plane  motions?  The  main 
assumption  that  needs  to  be  tested  is  whether  such  variations  produce  higher  order 
effects  as  we  assumed  in  this  thesis. 

2.  Is  there  a  significant  coupling  between  vertical  and  horizontal  plane  motions 
when  the  two  ships  are  under  towing?  In  order  to  answer  this  we  must  first 
produce  an  estimate  of  the  connection  forces  in  sway. 
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APPENDIX  A. 


REGULAR  WAVES  RESULTS 
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Figure  A-l.  Exciting  Forces  for  V=10kts,  1/L=0 . 1  and  (3=0° 


Figure  A-2.  Motions  of  the  Ships  for  V=10kts,  1/F=0. 1  and  |3=0° 
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Figure  A-3.  Exciting  Forces  for  V=10kts,  1/L=0. 1  and  |3=90° 
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Figure  A-7.  Force  Variation  for  Different  Headings  when  V=10kts  and  1/L=0. 1 
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Figure  A-9. 


KAIMALINO  Motion  for  Different  Headings  when  V=10kts  and 
1/L=0.1 


Figure  A- 10.  Connection  Force  for  Different  Towing  Lengths  when  V=10kts 

and  P=45° 
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Figure  A-l  1 .  SLICE  Motion  for  Different  Towing  Lengths  when  V=10kts  and 

(3=45° 


Figure  A-12.  KAIMALINO  Motion  for  Different  Towing  Lengths  when 

V=10kts  and  (3=45° 
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Figure  A- 13.  Connection  Forces  for  Different  Ship  Speeds  when  1/L=0. 1  and 

(3=45° 


T  [sec] 

Figure  A- 14.  SLICE  Motion  for  Different  Ship  Speeds  when  FL=0. 1  and  |3=450 

26 


4.5 
4 

3.5 
3 

2.5 

<£ 

2 

1.5 
1 

0.5 
0 

0  2  4  6  8  10  12  14 

T  [sec] 

Figure  A- 15.  KAIMALINO  Motion  for  Different  Ship  Speeds  when  1/L=0. 1  and 

(3=45° 


1 

1 

1 

1 

i 

i 

i 

-  -  5  w/o  con 
-  5  w/  con 

-  -  10  w/o  con 

-  10  w/  con 

15  w/o  con 
-  15  w /  con 

-  -  20  w/o  con 

if . 

_ A _ 

i 

i 

i 

i 

) 

i 

i 

i 

i 

r 

1 

_ 

1 

1 

i 

i 

i 

i 

-  2C 

w/  con 

1- 

) 

_ /. 

1 

1 

r 

( 

X 

s  I 

\l 

l 

l 

'  1  % 

!  \  \ 

1  A 
//  \ 

T 

i 

/ 

/  i 

i 

-  \Jr~n 

T  /V  n 

\t^^prz  ~  .  -  - i 

■ 

7  l  i  -  — f  i"  i  _ 

//  \  1  \  /  1  _ .  -  ~~vr  " 

T  Ws.'v.  J  >  '  ' 

vf  O  1 

_ 1 _ _ _ 1 _ 1 _ 

—  -  —  -  . 

27 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


28 


APPENDIX  B. 


RANDOM  WAVES  RESULTS 
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RMS  -n  RMS 


Heading 


igure  B- 1 . 


RMS  of  Connection  Force  vs.  Fleading  when  V=10kts  and 
1/L=0.01 


Figure  B-2.  RMS  of  Heave  Force  vs.  Heading  when  V=10kts  and  FL=0.01 
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Figure  B-3.  RMS  of  Vertical  Motion  vs.  Heading  when  V=10kts  and  1/L=0.01 


Figure  B-4.  Normalized  Force  (f/F3)  vs.  Heading  when  V=10kts  and  1/L=0.01 
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Heading 

Figure  B-5.  Vertical  Motion  vs.  Heading  when  V=10kts  and  1/L=0.01 


Heading 


Figure  B-6.  RMS  of  Connection  Force  vs.  Heading  when  V=10kts  and  1/L=0. 1 
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Figure  B-7.  RMS  of  Heave  Force  vs.  Heading  when  V=10kts  and  1/L=0.1 


Figure  B-8.  RMS  Value  of  Vertical  Motion  vs.  Heading  when  V=10kts  and 

1/LO.l 
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Heading 


Figure  B-9.  Normalized  Force  (f/F3)  vs.  Heading  as  a  Function  of  Sea  State 
when  V=10kts  and  1/L=0.1 


Heading 

Figure  B-10.  Vertical  Motion  vs.  Heading  when  V=10kts  and  1/L=0. 1 
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RMS  £  RMS 
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-11.  RMS  of  Connection  Force  vs.  Heading  as  a  Function  of  Sea  State 
when  V=10kts  and  1/L=1 


Figure  B-12.  RMS  of  Heave  Force  vs.  Heading  when  V=10kts  and  1/L=1 
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Figure  B-14.  Normalized  Force  (f/F3)  vs.  Heading  when  V=10kts  and  1/L=1 
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Figure  B-15.  Normalized  Vertical  Motion  vs.  Heading  when  V=10kts  and  1/L=1 


Figure  B-16.  RMS  of  Connection  Force  vs.  Heading  as  a  Function  of  1/L  when 

V=10kts 
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Figure  B-17.  RMS  of  Heave  Force  vs.  Heading  as  a  Function  of  1/L  when  V= 

lOkts 


Figure  B-18.  RMS  Value  of  Vertical  Motion  vs.  Heading  as  a  Function  of  1/L 

when  V=  lOkts 
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Figure  B-19.  Normalized  Force  (f/F3)  vs.  Fleading  as  a  Function  of  1/L  when 

V=10kts 


Figure  B-20.  Vertical  Motion  vs.  Heading  as  a  Function  of  1/L  when  V=10kts 
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Figure  B-2 1 .  RMS  of  Connection  Force  vs.  Heading  when  V=5kts  and  1/L=0. 1 


Figure  B-22.  RMS  of  Heave  Force  vs.  Heading  when  V=5kts  and  FL=0. 1 
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Heading 

Figure  B-23.  RMS  Value  of  Vertical  Motion  vs.  Heading  when  V=5kts  and 

1/L=0.1 


41 


Heading 


Figure  B-25.  Vertical  Motion  vs.  Heading  when  V=5kts  and  1/L=0. 1 


Figure  B-26.  RMS  of  Connection  Force  vs.  Heading  when  V=15kts  and  1/L=0. 1 
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Figure  B-28.  RMS  Value  of  Vertical  Motion  vs.  Heading  when  V=15kts  and 

1/LO.l 
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Figure  B-29.  Normalized  Force  (f/F3)  vs.  Heading  when  V=15kts  and  1/L=0. 1 


Figure  B-30.  Normalized  Vertical  Motion  vs.  Heading  when  V=15kts  and 

1/LO.l 
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Figure  B-3 1 .  RMS  of  Connection  Force  vs.  Heading  when  V=20kts  and  1/L=0. 1 


Figure  B-32.  RMS  of  Heave  Force  vs.  Heading  when  V=20kts  when  1/L=0. 1 
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Figure  B-33.  RMS  Value  of  Vertical  motion  vs.  Fleading  when  V=20kts  and 

1/L=0.1 
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Figure  B-35.  Vertical  Motion  vs.  Heading  when  V=20kts  and  1/L=0. 1 


Heading 


Figure  B-36.  RMS  of  Connection  Force  vs.  Heading  when  1/L=0. 1  and  Sea  State 

4 
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Figure  B-37.  RMS  of  Heave  Force  vs.  Heading  when  1/L=0. 1  and  Sea  State  4 


Heading 

Figure  B-38.  RMS  Value  of  Vertical  Motion  vs.  Heading  when  1/L=0. 1  and  Sea 

State  4 
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Figure  B-39.  Normalized  Force  (f/F3)  vs.  Heading  when  1/L=0.1  and  Sea  State  4 


Heading 

Figure  B-40.  Vertical  Motion  vs.  Heading  when  1/L=0. 1  and  Sea  State  4 
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APPENDIX  C. 


SEAKEEPING  RESULTS 
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Figure  C-l.  Sea  state  variation  when  V=10kts,  1/L=0 . 0 1  and  Sea  State  1 


Figure  C-2.  Sea  state  degradation  when  V=10kts,  1/L=0.0 1  and  Sea  State  1 
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Figure  C-3.  Sea  state  variation  when  V=10kts,  1/L=0.01  and  Sea  State  2 


Figure  C-4.  Sea  state  degradation  when  V=10kts,  1/L=0.0 1  and  Sea  State  2 
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Figure  C-5.  Sea  state  variation  when  V=10kts,  1/L=0.01  and  Sea  State  3 


Figure  C-6.  Sea  state  degradation  when  V=10kts,  1/L=0.0 1  and  Sea  State  3 


54 


90 


8 


Figure  C-7.  Sea  state  variation  when  V=10kts,  1/L=0 . 0 1  and  Sea  State  4 
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Figure  C-8.  Sea  state  degradation  when  V=10kts,  1/L=0.0 1  and  Sea  State  4 
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Figure  C-9.  Sea  state  variation  when  V=10kts,  1/L=0 . 0 1  and  Sea  State  5 


Figure  C-10.  Sea  state  degradation  when  V=10kts,  1/L=0.0 1  and  Sea  State  5 
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Figure  C-ll. 


Figure  C-12. 


Sea  state  variation  when  V=10kts,  1/L=0.01  and  Sea  State  6 
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Sea  state  degradation  when  V=10kts,  1/L=0.0 1  and  Sea  State  6 
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Figure  C-13. 


Figure  C-14. 


Sea  state  variation  when  V=10kts,  1/L=0.01  and  Sea  State  7 
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Sea  state  degradation  when  V=10kts,  1/L=0.0 1  and  Sea  State  7 


58 


90 


2.5 


Figure  C-15.  Sea  state  variation  when  V=10kts,  1/L=0.1  and  Sea  State  1 


Figure  C-16.  Sea  state  degradation  when  V=10kts,  1/L=0.1  and  Sea  State  1 
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Figure  C-17.  Sea  state  variation  when  V=10kts,  1/L=0.1  and  Sea  State  2 


Figure  C-18.  Sea  state  degradation  when  V=10kts,  1/L=0.1  and  Sea  State  2 
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Figure  C-19.  Sea  state  variation  when  V=10kts,  1/L=0.1  and  Sea  State  3 


Figure  C-20.  Sea  state  degradation  when  V=10kts,  1/L=0.1  and  Sea  State  3 
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Figure  C-2 1 .  Sea  state  variation  when  V=10kts,  1/L=0. 1  and  Sea  State  4 
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Figure  C-22.  Sea  state  degradation  when  V=10kts,  1/L=0.1  and  Sea  State  4 
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Figure  C-23.  Sea  state  variation  when  V=10kts,  1/L=0.1  and  Sea  State  5 


Figure  C-24.  Sea  state  degradation  when  V=10kts,  1/L=0.1  and  Sea  State  5 
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Figure  C-25.  Sea  state  variation  when  V=10kts,  1/L=0.1  and  Sea  State  6 


Figure  C-26.  Sea  state  degradation  when  V=10kts,  1/L=0.1  and  Sea  State  6 
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Figure  C-27.  Sea  state  variation  when  V=10kts,  1/L=0.1  and  Sea  State  7 
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Figure  C-28.  Sea  state  degradation  when  V=10kts,  1/L=0.1  and  Sea  State  7 
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Figure  C-29.  Sea  state  variation  when  V=10kts,  1/L=0.5  and  Sea  State  1 


Figure  C-30.  Sea  state  degradation  when  V=10kts,  1/L=0.5  and  Sea  State  1 
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Figure  C-3 1 .  Sea  state  variation  when  V=10kts,  1/L=0.5  and  Sea  State  2 


Figure  C-32.  Sea  state  degradation  when  V=10kts,  1/L=0.5  and  Sea  State  2 
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Figure  C-33.  Sea  state  variation  when  V=10kts,  1/L=0.5  and  Sea  State  3 


Figure  C-34.  Sea  state  degradation  when  V=10kts,  1/L=0.5  and  Sea  State  3 
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Figure  C-35.  Sea  state  variation  when  V=10kts,  1/L=0.5  and  Sea  State  4 
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Figure  C-36.  Sea  state  degradation  when  V=10kts,  1/L=0.5  and  Sea  State  4 
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Figure  C-37.  Sea  state  variation  when  V=10kts,  1/L=0.5  and  Sea  State  5 


Figure  C-38.  Sea  state  degradation  when  V=10kts,  1/L=0.5  and  Sea  State  5 
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Figure  C-39.  Sea  state  variation  when  V=10kts,  1/L=0.5  and  Sea  State  6 


Figure  C-40.  Sea  state  degradation  when  V=10kts,  1/L=0.5  and  Sea  State  6 
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Figure  C-41.  Sea  state  variation  when  V=10kts,  1/L=0.5  and  Sea  State  7 
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Figure  C-42.  Sea  state  degradation  when  V=10kts,  1/L=0.5  and  Sea  State  7 
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Figure  C-43.  Sea  state  variation  when  V=10kts,  1/L=1  and  Sea  State  1 


Figure  C-44.  Sea  state  degradation  when  V=10kts,  1/L=1  and  Sea  State  1 
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Figure  C-45.  Sea  state  variation  when  V=10kts,  1/L=1  and  Sea  State  2 


Figure  C-46.  Sea  state  degradation  when  V=10kts,  1/L=1  and  Sea  State  2 
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Figure  C-47.  Sea  state  variation  when  V=10kts,  1/L=1  and  Sea  State  3 


Figure  C-48.  Sea  state  degradation  when  V=10kts,  1/L=1  and  Sea  State  3 
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Figure  C-49.  Sea  state  variation  when  V=10kts,  1/L=1  and  Sea  State  4 
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Figure  C-50.  Sea  state  degradation  when  V=10kts,  1/L=1  and  Sea  State  4 
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Figure  C-5 1 .  Sea  state  variation  when  V=10kts,  1/L=1  and  Sea  State  5 


Figure  C-52.  Sea  state  degradation  when  V=10kts,  1/L=1  and  Sea  State  5 
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Figure  C-53.  Sea  state  variation  when  V=10kts,  1/L=1  and  Sea  State  6 


Figure  C-54.  Sea  state  degradation  when  V=10kts,  1/L=1  and  Sea  State  6 
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Figure  C-55.  Sea  state  variation  when  V=10kts,  1/L=1  and  Sea  State  7 
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Figure  C-56.  Sea  state  degradation  when  V=10kts,  1/L=1  and  Sea  State  7 
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APPENDIX  D.  HYDRODYNAMICS  DATA 


“MATDATA”  output  files  are  generated  by  SHIPMO  (Nash).  These  files  provide 
regular  wave  response  mass  added  and  excitation  force  matrix  constants  for  given  ship, 
speed,  and  wave  angle.  File  names  are  described  in  the  following  format: 

m  -matdata. 

s  or  k  -Vessel  simulated,  s-SLICE;  k-KAIMALINO. 
v  or  vh  -Motion  simulated,  v-vertical;  vh-all  six  degrees  of  freedom. 

Speed  -Zero  to  twenty  knots  in  one-knot  increments. 

Angle  -Zero  to  180  degrees  in  five-degree  increments. 


Example: 

rnkvh5_180.txt  =  Kaimalino,  motion  in  6-dof,  at  5  knots,  180°  wave  angle. 


SLICE  matdata  files 

KAIMALINO  matdata  files 

mkvhO  0.txt 

msvhO  0.txt 

mkvhl  0.txt 

msvhl  0.txt 

mkvh2  0.txt 

msvh2  0.txt 

mkvh3  0.txt 

msvh3  0.txt 

mkvh4  0.txt 

msvh4  0.txt 

mkvh5  0.txt 

msvh5  0.txt 

mkvh20  0.txt 

msvh20  0.txt 

mkvhO  5.txt 

msvhO  5.txt 

mkvhl  5.txt 

msvhl  5.txt 

mkvh2  5.txt 

msvh2  5.txt 

mkvh20  180.txt 

msvh20  180.txt 

(2  ships)  X  (21  speeds)  X  (37  angles)  =  1554  files 
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APPENDIX  E.  MATLAB  CODES  FOR  REGULAR  WAVES 


%  Vertical  Plane 

% 

clear 

% 

%  Get  run  info 

% 

V  =input(’Speed  (knots)  =  '); 
beta=input('Heading  (deg)  = '); 

1  =input('Length  (1/L)  =  '); 

% 

%  Get  tension  from  curvefitting  data. 

%  Applicable  for  speeds  between  1  and  20  ft/sec. 

% 

T=-1.762*VA4+63.675*VA3-580.8*VA2+2485.9*V-34.047; 

% 

V_string  =num2str(V); 
beta_string=num2str(beta); 

% 

%  The  matdata  output  files  default  to  the  vertical  only  format  when  the 
%  heading  angle  is  0  or  180  degrees. 

%  Set  up  file  reading  format. 

% 

trigg  =  30; 

f31oc  =  27;  f51oc=29; 
if  beta==0 
trigg  =  27; 

f31oc  =  26;  f51oc=27; 
elseif  beta==180 
trigg  =  27; 

f31oc  =  26;  f51oc=27; 
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end 

% 

%  Load  FRONT  SHIP  data  file  msvhV_beta.txt 

% 

load_filename=strcat(’msvh’,V_string,'_',beta_string,'.txt'); 

filename_s=load(load_filename); 

% 

%  Load  REAR  SHIP  data  file 

% 

load_filename=strcat('mkvh',V_string,'_',beta_string,'.txt'); 

filename_k=load(load_filename); 

% 

%  GENERAL  DATA 

% 

V=V*  1 .6878;  %  Convert  to  ft/sec 

lambda_min=20;  %  Min  wave  length  (ft) 

lambda_max=1000;  %  Max  wave  length  (ft) 

delta_lambda=20;  %  Wave  length  increment  (ft) 

rho=1.9905;  %  Water  density 

zeta=l ;  %  Regular  wave  height 

L=105;  %  Reference  length  for  nondimensionalization 

g=32.2;  %  Gravitational  constant 

x_s=-46;  %  FRONT  SHIP  attachment  point 

x_k=+40;  %  REAR  SHIP  attachment  point 

beta  =  beta*pi/180; 

lambda  =lambda_min:delta_lambda:lambda_max;  %  Vector  of  wavelengths 


wavenumber  =  2. 0*pi. /lambda; 

%  Wave  number 

omega 

=  sqrt(wavenumber*g); 

%  Wave  frequency 

omegae 

=  omega-wavenumber*V*cos(beta); 

%  Frequency  of  encounter 

period 

=  2.0*pi./omega; 

periode 

=  2. 0*pi. /omegae; 

omega 

=  omega'; 

omegae 

=  omegae'; 
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filesize  =  size(lambda); 
lambda_size=  trigg*filesize(2); 

% 

%  FRONT  SHIP 

% 

%  Set  mass  matrix  elements 
% 

M33s=filename_s(3:trigg:lambda_size,3); 

M35s=filename_s(3:trigg:lambda_size,5); 

M53s=filename_s(5:trigg:lambda_size,3); 

M55s=fdename_s(5:trigg:lambda_size,5); 

% 

%  Added  mass  terms 

% 

A3  3  s=filename_s(9 :  trigg :  lambda_size,3) ; 
A3  5  s=fdename_s(9 :  trigg :  lambda_size,5); 
A53s=fdename_s(l  1  :trigg:lambda_size,3); 
A55s=fdename_s(l  1  :trigg:lambda_size,5); 

% 

%  Damping  terms 

% 

B33s=fdename_s(  1 5  :trigg:lambda_size,3); 
B35s=filename_s(  1 5  :trigg:lambda_size,5); 
B53s=filename_s(17:trigg:lambda_size,3); 
B55s=fdename_s(17:trigg:lambda_size,5); 
% 

%  Hydrostatic  terms 
% 

C33s=filename_s(2 1  :trigg:lambda_size,3); 
C35s=filename_s(2 1  :trigg:lambda_size,5); 
C53s=fdename_s(23:trigg:lambda_size,3); 
C55s=fdename_s(23:trigg:lambda_size,5); 
% 
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%  Total  exciting  forces 

% 

F3s_t_amp=filename_s(f31oc:trigg:lambda_size,5); 
F5  s_t_amp=filename_s(f5  loc :  trigg :  lambda_size,5); 
F3  s_t_pha=fdename_s(f3  loc :  trigg :  lambda_size  ,6) ; 
F  5  s_t_pha=fdename_s(f5  loc :  trigg :  lambdasize  ,6) ; 
F3s_t=F3s_t_amp.*exp(i*F3s_t_pha.*pi/180.0); 

F5  s_t=F5  s_t_amp .  *  exp(i  *  F5  s_t_pha .  *  pi/ 1 8  0 . 0 ) ; 

% 

%  Froude/Krylov  exciting  forces 
% 

F3s_f_amp=filename_s(f31oc:trigg:lambda_size,l); 
F  5  s_f_amp=filename_s(f5  loc :  trigg:  lambdasize,  1 ); 
F3  s_f_pha=fdename_s(f3  loc :  trigg  :lambda_size,2) ; 
F  5  s_f_pha=fdename_s(f5  loc :  trigg  :lambda_size,2) ; 
F3s_f=F3s_f_amp.*exp(i*F3s_f_pha.*pi/180.0); 
F5s_f=F5s_f_amp.*exp(i*F5s_f_pha.*pi/180.0); 

% 

%  Diffraction  exciting  forces 
% 

F3s_d_amp=filename_s(f31oc:trigg:lambda_size,3); 
F5s_d_amp=filename_s(f51oc:trigg:lambda_size,3); 
F3s_d_pha=fdename_s(f31oc:trigg:lambda_size,4); 
F  5  s_d_pha=fdename_s(f5  loc :  trigg :  lambda_size,4); 
F3s_d=F3s_d_amp.*exp(iH5F3s_d_pha.*pi/180.0); 
F5s_d=F5s_d_amp.*exp(i*F5s_d_pha.*pi/180.0); 
% 

%  REAR  SHIP 

% 

%  Set  mass  matrix  elements 
% 

M33k=fdename_k(3:trigg:lambda_size,3); 

M35k=filename_k(3:trigg:lambda_size,5); 
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M53k=filename_k(5:trigg:lambda_size,3); 

M55k=filename_k(5:trigg:lambda_size,5); 

% 

%  Added  mass  terms 

% 

A33k=filename_k(9:trigg:lambda_size,3); 
A35k=filename_k(9:trigg:lambda_size,5); 
A53k=filename_k(l  1  :trigg:lambda_size,3); 
A55k=filename_k(  1 1  :trigg:lambda_size,5); 

% 

%  Damping  terms 

% 

B  3  3k=filename_k(  1 5 :  trigg :  lambdasize,  3  ) ; 
B35k=filename_k(  1 5  :trigg:lambda_size,5); 

B  5  3k=filename_k(  1 7 :  trigg :  lambdasize,  3  ) ; 
B55k=filename_k(17:trigg:lambda_size,5); 

% 

%  Hydrostatic  terms 

% 

C33k=filename_k(2 1  :trigg:lambda_size,3); 
C35k=niename_k(2 1  :trigg:lambda_size,5); 
C53k=fdename_k(23:trigg:lambda_size,3); 

C  5  5k=fdename_k(2  3 :  trigg :  lambdasize,  5  ) ; 

% 

%  Total  exciting  forces 

% 

F3k_t_amp=filename_k(f31oc:trigg:lambda_size,5); 
F5k_t_amp=filename_k(f51oc:trigg:lambda_size,5); 
F3  k_t_pha=n  lename_k((3  loc:  trigg:  lambda_sizc,6); 
F5k_t_pha=filename_k(f5  loc :  trigg:  lambda_size,6) ; 
F3k_t=F3k_t_amp.*exp(i*F3k_t_pha.*pi/180.0); 
F5k_t=F5k_t_amp.H5exp(i*F5k_t_pha.*pi/180.0); 

% 
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%  Froude/Krylov  exciting  forces 
% 

F3k_f_amp=filename_k(f3  loc :  trigg :  lamb  dasize ,  1 ) ; 

F5k_f_amp=filename_k(f5  loc :  trigg :  lamb  dasize ,  1 ) ; 
F3k_f_pha=filename_k(f31oc:trigg:lambda_size,2); 
F5k_f_pha=filename_k(f51oc:trigg:lambda_size,2); 
F3k_f=F3k_f_amp.*exp(i*F3k_f_pha.*pi/180.0); 

F  5  k_f=  F  5k_f_amp .  *  exp(i  *  F5k_  f_p  h  a .  *  p  i  / 1 8  0 . 0 ) ; 

% 

%  Diffraction  exciting  forces 
% 

F3k_d_amp=filename_k(01oc:trigg:lambda_size,3); 

F5k_d_amp=filename_k(f51oc:trigg:lambda_size,3); 

F3k_d_pha=filename_k(f3  loc :  trigg :  lambda_size,4); 

F5k_d_pha=fdename_k(f5  loc :  trigg :  lambda_size,4); 

F3k_d=F3k_d_amp.*exp(i*F3k_d_pha.*pi/180.0); 

F5k_d=F5k_d_amp.*exp(i*F5k_d_pha.*pi/180.0); 

% 

%  MATCHING  CONDITION 

% 

A33bar_s=-(omegae.A2).*(M33s+A33s)+i*omegae.*B33s+C33s; 

A35bar_s=-(omegae.A2).*(M35s+A35s)+i*omegae.*B35s+C35s; 

A53bar_s=-(omegae.A2).*(M53s+A53s)+i*omegae.*B53s+C53s; 

A55bar_s=-(omegae.A2).*(M55s+A55s)+i*omegae.*B55s+C55s; 

A33bar_k=-(omegae.A2).*(M33k+A33k)+i*omegae.*B33k+C33k; 

A35bar_k=-(omegae.A2).*(M35k+A35k)+i*omegae.*B35k+C35k; 

A53bar_k=-(omegae.A2).*(M53k+A53k)+i*omegae.*B53k+C53k; 

A55bar_k=-(omegae.A2).H5(M55k+A55k)+i*omegae.*B55k+C55k; 

% 

mu3_s=(A55bar_s.*F3s_t-A35bar_s.*F5s_t)./(A33bar_s.*A55bar_s- 

A35bar_s.*A53bar_s); 

nu3_s=(A55bar_s+A35bar_s*x_s)./(A33bar_s.*A55bar_s-A35bar_s.*A53bar_s); 
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mu5_s=(A53bar_s.*F3s_t-A33bar_s.*F5s_t)./(A53bar_s.*A35bar_s- 

A33bar_s.*A55bar_s); 

nu5_s=(A53bar_s+A33bar_s*x_s)./(A53bar_s.*A35bar_s-A33bar_s.*A55bar_s); 

mu3_k=(A55bar_k.*F3k_t-A35bar_k.*F5k_t)./(A33bar_k.*A55bar_k- 
A3  5bar_k.  *  A53bar_k); 

nu3_k=(A55bar_k+A35bar_k*x_k)./(A33bar_k.*A55bar_k- 
A3  5bar_k.  *  A53bar_k); 

mu5_k=(A53bar_k.*F3k_t-A33bar_k.*F5k_t)./(A53bar_k.*A35bar_k- 

A33bar_k.*A55bar_k); 

nu5_k=(A53bar_k+A33bar_k*x_k)./(A53bar_k.*A35bar_k- 

A33bar_k.*A55bar_k); 

% 

a=mu3_s -mu5_s  *  x_s-mu3_k+mu5_k  *  x_k ; 
b=nu3_s-nu5_s*x_s+nu3_k-nu5  k*x  k; 
f=a./(l/T+b); 

% 

f_s=-f; 
f_k=f; 

eta3_s=mu3_s+nu3_s.  *f_s; 
eta5_s=mu5_s+nu5_s.  *f_s; 
eta3_k=mu3_k+nu3_k.  *  f_k; 
eta5_k=mu5_k+nu5_k.  *f_k; 
xi_s=eta3_s-eta5_s*x_s; 
xi_k=eta3_k-eta5_k*x_k; 
xiO_s=mu3_s-mu5_s*x_s; 
xi0_k=mu3_k-mu5_k*x_k; 

% 

%  SAMPLE  FIGURES 

% 

figure  (1) 

plot(period,abs(xiO_s),'r’,  period,  abs(xi_s),’b'),  grid,  legend('w/o  connection’,'with 

connection’) 

title('Leading  Ship  Motion’), xlabel(’T  [sec]’),ylabel('\xi_V  [ft/ft]') 


%  Connection  force  on  FRONT  SHIP 
%  Connection  force  on  REAR  SHIP 
%  FRONT  SHIP  heave 
%  FRONT  SHIP  pitch 
%  REAR  SHIP  heave 
%  REAR  SHIP  pitch 
%  FRONT  SHIP  motion  at  connection 
%  REAR  SHIP  motion  at  connection 
%  FRONT  SHIP  motion  at  connection  for  zero  f 
%  REAR  SHIP  motion  at  connection  for  zero  f 
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figure  (2) 

plot(period,abs(xiO_k),'r’,period,abs(xi_k),’b'),grid,legend('w/o  connection',' with 
connection') 

title('Trailing  Ship  Motion’), xlabel('T  [sec]’),ylabel('\xi_V  [ft/ft]’) 


figure  (3) 

plot(period,abs(f),'r’,  period,  F3s_t_amp,’b’,  period,  F3k_t_amp,’g'),  grid,  legend('conn 
ection  force', 'wave,  leading  ship’, 'wave,  trailing  ship’) 

title('Exciting  Forces'), xlabel(’T  [sec]’),ylabel('F  [lbs/ft]’) 
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APPENDIX  F.  MATLAB  CODES  FOR  RANDOM  WAVES 


%  This  section  is  added  to  the  codes  for  regular  waves 
hs=[.85  2.25  4.25  6.75  10  17  32.5]; 

% 

%  Begin  random  wave  calculations 
%  Fully  developed  seas  -  PM  spectrum 

% 

iHS=0; 

for  1=1 : 7;  %  Loop  on  significant  wave  height 

HS=hs(I); 
iHS=iHS+l; 

% 

%  Random  wave  calculations 
%  Pierson-Moscowitz  spectrum 

% 

POWER  =-.032*(g/HS)A2; 

S  =(0.008  l*gA2).*exp(POWER./(omega.A4))./(omega.A5); 

Se  =S./(l-(2.0/g)*omega*V*cos(beta));  %  Convert  S(w)  to  S(we) 

% 

%  Define  response  spectra 
% 

Sf  =((abs(f)).A2).*Se; 

Sxi_s  =((abs(xi_s)).A2).*Se; 

Sxi_k  =((abs(xi_k)).A2).*Se; 

Sxi0_s  =((abs(xiO_s)).A2).*Se; 

Sxi0_k  =((abs(xiO_k)).A2).*Se; 

SF3s_t  =((abs(F3s_t)).A2).*Se; 

SF3k_t  =((abs(F3k_t)).A2).*Se; 

% 

%  Initializations 

% 

Sf_i=0; 
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Sxi_s_i=0; 

Sxi_k_i=0; 

Sxi0_s_i=0; 

Sxi0_k_i=0; 

SF3s_t_i=0; 

SF3k_t_i=0; 

% 

%  Integral  S(w)*|RA0|A2 
% 


omegae(I))); 


for  1=2: 1  :filesize(2), 

Sf_i=Sf_i+0.5*abs((Sf(I)+Sf(I-l)))*abs((omegae(I- 1  )- 


Sxi_s_i=Sxi_s_I+0.5*abs((Sxi_s(I)+Sxi_s(I- 
1  )))*  abs((omegae(I- 1  )-omegae(I))); 

Sxi_k_i=Sxi_k_i+0.5*abs((Sxi_k(I)+Sxi_k(I- 

l)))*abs((omegae(I-l)-omegae(I))); 

SxiO_s_i=  SxiO_s_i  +  0.5*abs((Sxi0_s(I)  +  SxiO_s(I-l)))  * 
abs((omegae(I- 1  )-omegae(I))); 

SxiO_k_i=  SxiO_k_i  +  0.5*abs((Sxi0_k(I)  +  SxiO_k(I-l)))  * 
abs((omegae(I- 1  )-omegae(I))); 

SF3s_t_i=  SF3s_t_i  +  0.5*abs((SF3s_t(I)  +  SF3s_t(I-l )))  * 
abs((omegae(I- 1  )-omegae(I))); 

SF3k_t_i=  SF3k_t_i  +  0.5*abs((SF3k_t(I)  +  SF3k_t(I-l))) 
*  abs((omegae(I-l)-omegae(I))); 

end 


% 

%  RMS  values 
% 

RMS_f  =  sqrt(Sf_i); 
RMS_xi_s  =  sqrt(Sxi_s_i); 
RMS_xi_k  =  sqrt(Sxi_k_i); 
RMS_xiO_s  =  sqrt(SxiO_s_i); 
RMSxiOk  =  sqrt(SxiOki); 
RMS_F3s_t  =  sqrt(SF3s_t_i); 
RMS_F3k_t  =  sqrt(SF3k_t_i); 
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vHS(ibeta,iHS)=HS; 

RMS_fF(ibeta,iHS)=RMS_f/RMS_F3s_t; 

RMS_xX(ibeta,iHS)=RMS_xi_s/RMS_xiO_s; 

RMS_ff(ibeta,iHS)=RMS_f; 

RMS_xx(ibeta,iHS)=RMS_xi_s; 

RMS_FF3s_t(ibeta,iHS)=RMS_F3s_t; 
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APPENDIX  G.  MATLAB  CODES  FOR  SEAKEEPING 

EVALUATION 

%  This  section  is  added  at  the  end  the  codes  for  random  waves 
RMS_con=RMS_ff+RMS_FF3s_t; 

RMS_uncon=RMS_FF3s_t; 

%  This  section  is  a  sample  for  getting  sea  state  figures  of  Appendix.  C  after  you 
build  a  matrix  of  RMS  values  above  for  different  conditions 

t=linspace(0,2*pi,73)'; 

figure(l) 

polar(t,SS  1005(:, 23)), hold  on 
polar(t,7  *  ones(73 , 1 )) 
figure(2) 

polar(t,SS  10050, 22)), hold  on 
polar(t,6  *  ones(73 , 1 )) 
figure(3) 

polar(t,SS  10050,21)), hold  on 
polar(t,5  *  ones(73 , 1 )) 
figure(4) 

polar(t,SS1005(:,20)),hold  on 
polar(t,4  *  ones(73 , 1 )) 
figure(5) 

polar(t,SS1005(:,19)),hold  on 
polar(t,3  *  ones(73 , 1 )) 
figure(6) 

polar(t,SS  1 005(:,  1 8)), hold  on 
polar(t,2  *  ones(73 , 1 )) 
figure(7) 

polar(t,SS1005(:,17)),hold  on 
polar(t,  1  *  ones(73 , 1 )) 
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